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Abstract. This review consists of two parts. The first part establishes certain astrophysical bounds 
on the smoothness of classical spacetime. Some of the best bounds to date are based on the absence 
of vacuum Cherenkov radiation in ultrahigh-energy cosmic rays. The second part discusses possible 
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fundamental length scale of quantum spacetime may be different from the Planck length. 
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1. INTRODUCTION 

^ . The present contribution addresses the following basic question: does space remain 

smooth as one probes smaller and smaller distances ? A conservative limit on the typical 
length scale i of any nontrivial small-scale structure of space results from the fact that 
^ . the experimental data from particle accelerators such as LEP at CERN and the Tevatron 

t^. \ at Fermilab are perfectly well described by relativistic quantum fields over a smooth 

manifold (specifically, Minkowski spacetime), 

O i £ |LEP/Tevatron < 1Q _ 18 m ^ hcj{ ^^ (1) 

r- : 

Remark that the last approximate equality in (1) follows from the Heisenberg uncertainty 
principle applied to the matter probes, even though the bounded spacetime length £ itself 
may be a purely classical quantity. 

Yet, astrophysics provides us with very much higher energies to probe spacetime. A 
possible strategy is then 

• to consider the phenomenology of simple spacetime models; 

• to obtain bounds on the model parameters from ultrahigh-energy cosmic rays; 

• to establish the main theoretical implications. 

Some of our recent work has pursued the above strategy and the aim of the present 
article is to review this work in a coherent fashion. In Sec. 2, we discuss the phenomenol- 
ogy of two simple models for modified photon propagation [1, 2]. In Sees. 3.1 and 3.2, 
we obtain ultrahigh-energy cosmic-ray (UHECR) bounds [1, 2, 3] on the parameters of 
the two types of models considered and, in Sec. 3.3, we discuss some implications. In 
Sec. 4, which can be read independently of Sees. 2 and 3, we put forward two conjec- 
tures [4] on the fundamental length of a hypothetical small-scale structure of quantum 
spacetime. In Sec. 5, we summarize our results. 
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2. PHENOMENOLOGY 



Two simple models of photon propagation will be presented in this section. The first 
model has the standard quadratic Maxwell action density term integrated over a flat 
spacetime manifold with "defects." For this first model, the modified photon dispersion 
relation will be calculated in the long-wavelength approximation. 

The second model has a quadratic modified-Maxwell term integrated over standard 
Minkowski spacetime. The corresponding photon dispersion relation can be readily 
obtained. Having standard Dirac fermions coupled to photons with modified propaga- 
tion properties, the process of "vacuum Cherenkov radiation" may be allowed for cer- 
tain combinations of Lorentz- violating parameters. For the second model, the relevant 
Cherenkov energy threshold will be given explicitly. 



2.1. Model 1 

Consider standard Quantum Electrodynamics (QED) [5] over a classical spacetime- 
foam 1 manifold (with details to be specified later), which has the following action: 

S Mod e\i=f f ' d 4 x (-ir 7 ^r 7 VCT F MV WFp (7 W + v7W(/(i^- e A M W)-M)^)), 

(2) 

where F j _ lv (x) = d ll A v (x) — d v A^(x) is the Maxwell field strength tensor of the gauge 
field A^ \x) and yf(jc) the Dirac spinor field corresponding to, e.g., a proton with electric 
charge e and rest mass M. 

The prime on the integral sign of (2) indicates the restriction to the flat Minkowski 
part of the manifold with specific boundary conditions at certain submanifolds called 
"defects." The procedure can perhaps best be illustrated by the example of Fig. 1 on the 
next page, where the static "wormholes" [6, 7] are sliced off (or, more accurately, the 
wormhole "throats" are taken to have zero lengths) and the resulting holes/defects in 
flat spacetime are given appropriate boundary conditions for the vector and spinor fields 
appearing in (2); see Fig. 2 for a sketch. 

The coordinates of flat spacetime are denoted by (x^) = (x , x) = (ct, X . X , X 3 ) and 
the standard Minkowski metric is (r]^ v ) = diag(l,— 1,— 1,— 1). The Levi-Civita symbol 
£nvpa, which will be used later, is normalized by £0123 = 1. The direction of a 3-vector 
x is given by the unit 3-vector x = x/|x|. In this section and the next, we use natural 
units with c = fi= 1, but, occasionally, we display c or fi in order to clarify the physical 
dimensions of a particular expression. The dimensionful constants c and h will, however, 
occupy the center of the stage in Sec. 4. 



A note of caution may be in order, as the terminology "spacetime foam" is often considered to refer 
solely to the quantum structure of spacetime [6, 7, 8, 9], whereas, here, the picturesque designation 
"classical spacetime-foam manifold" simply refers to a classical manifold with nontrivial small-scale 
structure (resembling, for example, a well-known Swiss cheese). The elusive "quantum structure of 
spacetime" will be discussed further in Sec. 4, whereas we remain with the more or less familiar classical 
spacetime in Sees. 2 and 3. 
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FIGURE 1. Example of a static (time-independent) classical spacetime foam, with one spatial dimen- 
sion suppressed. The corresponding spacetime is mostly flat (Minkowski-like) and has several permanent 
wormholes added. 



Next, take some very simple classical spacetime-foam models with: 

• identical defects (spacetime size b) embedded in Minkowski spacetime, 

/ -4 \ 

• a homogeneous and isotropic distribution of defects (spacetime density n= l/l ) , 

• a strong dilution of defects (b<til), 

where the last requirement, in particular, is a purely technical assumption in order to 
simplify the calculation. 

Remark that, for the classical spacetime shown in Fig. 1, there is no topology change 
(in fact, it is unclear whether or not topology change is allowed at all; see, e.g., Ref. [9]). 
As far as our calculations are concerned, the detailed dynamics of the classical defects 
considered does not appear to be important; what matters are average quantities such as 
the typical defect size and separation. 

Now calculate the proton and photon dispersion relations in the long-wavelength 
approximation, X ^>1^> b. No details of the calculation [1] will be given here but only 
the heuristics for the photon case: localized defects correspond to fictional multipoles 
which affect the propagation of electromagnetic plane waves. Electromagnetic- wave 



FIGURE 2. Three-space from a single wormhole-like defect (two spheres with equal radius b and 
distance d >2b between their centers) embedded in K 3 , with the "interiors" of the two spheres removed 
and their points identified by reflection in a central plane with normal unit vector a. If the "long distance" 
d between the centers of the wormhole mouths is of the same order of magnitude as their width 2 b, the 
defect is effectively localized in the ambient three-space. 
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propagation over a spacetime with defects is then different from propagation over a 
perfectly smooth Minkowski spacetime. In fact, the basic physics is the same as that of 
the so-called "Bethe holes" in waveguides [10] . It is clear that different types of localized 
defects give essentially the same result, just with different numerical coefficients. 

The modified proton (p) and photon (7) dispersion relations for wave number k = 
|k| =2k/X can, therefore, be given in completely general form (that is, independent of 
the detailed calculations): 

CO 2 = m z p c 4 /h 2 + c 2 p k 2 + 0(c 2 p b 2 k 4 ) , (3a) 

co 2 = (I + 02F) c 2 p k 2 + a 4 F c 2 p b 2 k 4 + 0(c 2 p b 4 k 6 ), (3b) 

with c 2 simply defined as the coefficient of the quadratic proton term, effective on/off 

factors 02,04 e { — 1,0, +1}, an effective defect size b, and an effective excluded- 
volume factor _ _ 

F={b/l)\ (4) 

which, for the moment, is assumed to be much less than unity. 

Observe that the calculated dispersion relations as given by (3ab) do not contain cubic 
terms in k, consistent with general arguments based on coordinate independence and ro- 
tational invariance [11]. Furthermore, the photon dispersion relations found are the same 
for both polarization modes (absence of birefringence) because of the assumed isotropy 
of the defect distributions. However, the photon dispersion relations do show birefrin- 
gence, but still no cubic terms, if there is assumed to be an anisotropic distribution of 
intrinsically asymmetric defects (having, for example, defects from Fig. 2 with all axes 
a aligned). 

As mentioned above, calculations of certain simple spacetime-foam models give 
the effective dispersion-relation parameters (with tildes) in terms of the underlying 
spacetime parameters (with bars): 

b = pb, T=xh o 2 = -l, 04 = 1 , (5) 

where /3 and % are positive numerical constants of order unity, which depend on the 
details of the defects considered. Note that, for static defects as in Fig. 2, with size b 
and separation /, a physically more appropriate definition of F would be (b/l) 3 , but, 
mathematically, definition (4) with (b/l) 4 can still be used. 

With intrinsically different boundary conditions at the defect locations for the different 
types of matter fields (spinor and vector), the quadratic terms of the proton and photon 
dispersion relations can be expected to be different in general. Having defined the 
quadratic proton term in (3a) as c 2 k 2 , the quadratic photon term in (3b) will then differ 

from c 2 k 2 , unless the defects are infinitely small or infinitely far apart. The resulting 
unequal proton and photon velocities may lead to new types of decay processes [12], for 
example, vacuum Cherenkov radiation p — > p 7 for the case of a negative coefficient 02F 
in (3b). This particular decay process will be discussed further in the next subsection. 

Regardless of the origin and interpretation, it is clearly important to obtain bounds on 
the effective parameters F and b in the modified photon dispersion relation (3b) for a 
proton dispersion relation given by (3a). 
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2.2. Model 2 



In the previous subsection, we have found a quadratic photon term of the dispersion 
relation which was modified by a single parameter F determined by the underlying struc- 
ture of spacetime. This type of modified photon propagation can easily be generalized. 

Consider the following action for a Lorentz- violating deformation of QED: 

S|\/lodel2 = SmodM + S s tandD, (6) 

with a modified-Maxwell term [13, 14], 

S mod M = J r /x (- i {r\Wr\™ + K^P°)F^ v (x)F p(T (x)) , (7) 
and the standard Dirac term [5] for a spin-^ particle with charge e and mass M, 

SstandD = J d 4 * (f (i ^ - eAp (*)) - M) ¥ (x) . (8) 

Theory (6) is gauge-invariant, CPT-even, and power-counting renormalizable. 

The quantity K^ vpa in the modified-Maxwell term (7) is a constant background tensor 
with real and dimensionless components. This background tensor K^ vpa has, in fact, 
the same symmetries as the Riemann curvature tensor and a double trace condition 
K^^y — 0, so that there are 20 — 1 = 19 independent components. All components of the 
fc-tensor in (7) are assumed to be sufficiently small in order to ensure energy positivity. 

As the ten birefringent parameters are already constrained at the 10~ 32 level [15], 
restrict the theory to the nonbirefringent sector by making the following Ansatz [16]: 

for a symmetric and traceless matrix K^ v with 10—1 =9 independent components. 
Rewrite these nine Lorentz- violating "deformation parameters" K^ v as follows: 

(^ v )=diag(l,i,i,i)K 00 +(5K^), 5k 00 = 0, (10) 

with one independent parameter 7c 00 for the spatially isotropic part of K^ v and eight 
independent parameters 8k^ v which need not be smaller than 7c 00 . 
For later use, also define a vector a in parameter space R 9 : 
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FIGURE 3. Feynman diagram for vacuum Cherenkov radiation. 



which is taken to have the standard Euclidean norm, 

i«i 2 =i(«') 2 - 

1=0 

Note that the negative of the isotropic parameter a corresponds, in leading order, to the 
coefficient 02 F of the modified photon dispersion relation (3b) calculated for Model 1. 

It may be of interest to mention that, whereas Model 2 has constant deformation 
parameters (9), another model has a small stochastic parameter g(x) multiplying a CPT- 
odd action density term e^ vpa F^ v {x)F pa {x). This particular stochastic model has been 
studied in the long-wavelength approximation in Ref. [17]. Recently, we have become 
aware of a complementary paper [18], which considers similar models in the short- 
wavelength approximation. 

Model 2 has, for appropriate deformation parameters a, a maximum proton velocity 
larger than the phase velocity of light, which allows for vacuum Cherenkov radiation 
p — > py. This particular decay process has been studied classically by Altschul [19] 
and quantum-mechanically at tree-level (Fig. 3) by Kaufhold and the present author 
[2, 12]. The radiated-energy rate of a primary (on-shell) particle with point charge 
ZpHm e, mass M pr \ m > 0, momentum q pr i m , and ultrarelativistic energy E pr \ m ~ c |q pr 
is asymptotically given by [2] 



dWModel2(qprim, £prim) 

dt 



Zprim ^ % (Qprim) E prlm /h > (^) 



2 



c 'prim^ c thresh 



with a nonnegative dimensionless coefficient £ (q pr im) from appropriate contractions of 
the K-tensor with two rescaled g-vectors and the squared threshold energy [19] 



M 2 ■ c 4 

El u = + 0(M 2 - c 4 ) (14) 

^thresh „ / n , ; , ~ ,t x ^ ^^prim c )i \ L *) 

R £ (a° + aJ q J - + W k q J - q; „■ ) 

c V ^pnm ^pnm ^prim / 



in terms of a regularized ramp function R e (x) = £ + (x + \x\ )/2 for a real variable x and 
an arbitrarily small positive parameter £. 

Exact tree-level results have been obtained recently [20] for the restricted model 
(labeled "isotropic case") with (Xq = a > and a 1 = 0, for / = 1, . . . ,8. Setting again 
h = c = 1, the radiated-energy rate for a spin-^ Dirac particle (charge e, mass M, and 
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energy E above threshold) is given by the following expression: 



.^isotropiccase 2 



dt An 3a~;EVE 2 -M 2 

x|2(a^ + 4a + 6) E 2 - (2 + ceo) 

x ^3{l + ao)M 2 + 2{3 + 2ao)^^^EVE^I^y (15) 
The high-energy expansion of (15) for fixed parameters 0Cq and M reads 

iTrrisotropiccase 2 
flKV Model2 



dt 



+ (- 1 + ^ a »-| a » 2 + 0(o » 3) ) M2 + (^ 1 >' (16) 

which shows the quadratic behavior of (13) with factor Z 2 r - Im = 1 and constant coefficient 
t, = (7/24) ceo. From (15), one also obtains the exact threshold energy (temporarily 
reinstating c): 

r^Model 2. isotropic case Mc r jz /n\ 

thresh =-= V / lW2, (17) 



which reproduces (14) for M — M prim and small enough positive ceo = 05°. Incidentally, 
the numerical value of ceo cannot be too large, as the factors \J2. — ceo in (15) make clear. 

Returning to Model 2 with generic deformation parameters ce z , it is to be expected that 
vacuum Cherenkov radiation only occurs for those parameters a for which the phase 
velocity of light is less than the maximal attainable velocity c of the charged particle in 
theory (6), v p h < c. In fact, this phase- velocity condition corresponds to having a positive 
argument of the ramp function on the right-hand side of (14). The relevant domain in 
parameter space is given by 

d{ Z^ = {^^ 9 - v Xe«3 (a° + a J # + fi'"W) >0}, (18) 

for arbitrary unit 3-vector x = x/|x|. The superscript '(open)' in (18) refers to the use 
of the open relation symbol '>' on the right-hand side instead of the closed symbol '>', 
because there is no vacuum Cherenkov radiation for the case of v p h = c. Most likely, the 
domain (18) constitutes a significant part of the physical domain of theory (6), where, 
e.g., microcausality holds. Hence, the subscript 'causal' in (18). 

Microcausality of Lorentz-violating theories has been discussed in, for example, 
Refs. [21, 22, 23]. To this author, it is clear that violation of microcausality is unac- 
ceptable in the type of theories considered here, because, for a spacelike separation 
of two events, the time order can be interchanged by an appropriate observer Lorentz 
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transformation. Perhaps not unacceptable (at least, according to our current experimen- 
tal knowledge) is a new class of particle instabilities which occur at ultralarge three- 
momentum due to energy non-positivity (at least, in "concordant" frames, where the 
Lorentz- violating parameters are relatively small). 

Let us end this subsection with a general remark on Model 2 (see also Sec. 3.3 for 
further discussion). New phenomena at the energy scale .Epianck ~ 10 19 GeV [6, 7] may 
lead to Lorentz violation in the low-energy theory, partially described by the model 
action (6). But the resulting Lorentz- violating parameters a need not be extremely 
small (e.g., suppressed by inverse powers of ^pianck) an d can be of order unity, as 
long as the theory remains physically consistent. In fact, the modified-dispersion-relation 
calculations from the previous subsection provide an example: the quadratic coefficient 
F in (3b) can, in principle, be close to one. Hence, it is important to obtain as strong 
bounds as possible on all deformation parameters a. 

3. UHECR CHERENKOV BOUNDS 

The goal of this section is to establish bounds on the parameters of the two Lorentz- 
violating models discussed in the previous section. One type of bounds relies on the 
process of vacuum Cherenkov radiation (Fig. 3) already mentioned in Sec. 2.2. 
The basic idea is simple [24, 25]: 

• if vacuum Cherenkov radiation has a threshold energy thresh (b,l, cc) and the ra- 
diation rate above threshold is not suppressed, then UHECRs with E prim > thresh 
cannot travel far (certainly not over distances of the order of megaparsecs), as they 
rapidly radiate away their energy; 

• observing an UHECR then implies a primary energy E pr \ m at or below threshold, 
which gives bounds on combinations of b, I, and a. 

Expanding on the last point, the following inequality must hold for the measured primary 
energy E pr \ m of an UHECR compared to the theoretical result for the threshold energy: 

£prim <£thresh(>,^ «), (19) 

which can be written as an upper bound on the parameters of the Lorentz-violating 
theory considered (here, taken to be b, I, and a, from the two models discussed in the 
two previous subsections). 

Incidentally, other types of bounds [1] have been obtained (specifically, from the lack 
of time dispersion and Rayleigh-like scattering in a particular TeV gamma-ray flare from 
the active galaxy Mkn 421), but the resulting bounds are, for the moment, less tight than 
those from the inferred absence of vacuum Cherenkov radiation. 



3.1. UHECR Bounds - Model 1 

In order to obtain Cherenkov bounds [1] for the two parameters b and I of Model 1, 
start with the spectacular event [26] shown in Fig. 4 a few pages later. This particular 
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Extensive Air Shower (EAS) event corresponds, in fact, to the most energetic particle 
observed to date, having an energy E pr \ m ss 300 EeV = 3 x 10 11 GeV = 3 x 10 20 eV. 

From the EAS event of Fig. 4, we obtain the following Cherenkov-like bounds on the 
quadratic and quartic coefficients of dispersion relation (3b): 

-3xl0- 23 < d 2 F < 3xl0- 23 , (20a) 

-(7xl0~ 39 m) 2 < o 4 Fb 2 < (5xl0" 38 m) 2 , (20b) 

based on the detailed analysis of Ref. [27] (some back-of-the-envelope calculations have 
been presented in App. B of Ref. [28]). For bounds (20ab), the primary was assumed to 
be a proton (p) with standard partonic distributions and, as mentioned before, its energy 
was determined to be E p ks 3 x 10 11 GeV. 

Note that the above bounds are two-sided, whereas genuine Cherenkov radiation 
would only give one-sided bounds (the electromagnetic-wave phase velocity must be 
less than the maximum particle velocity). The reason for having two-sided bounds is 
that, in addition to Cherenkov radiation, another type of process can occur, namely, 
proton break-up p — > pe + e~ (pair-production by a virtual gauge boson), which gives 
the "other" sides of the bounds [27]. 

It is important to understand the dependence of bounds (20ab) on the assumed primary 
mass (here, taken as M pr j m = 0.94 GeV/c 2 ) and primary energy (here, taken as E pr \ m = 
3 x 10 11 GeV). The high and low limits of bounds (20a) and (20b) are multiplied by 
scaling factors f a high/low an d fb high/low respectively, which are approximately given 

by 

/a,hi g h ~ /a,bw » (M prim c 2 /0.94GeV) 2 (3 x 10 11 GeV/£ prim ) 2 , (21a) 
/b,hi g h ~ /b,iow ~ (M prim c 2 /0.94GeV) 2 (3 x 10 11 GeV/£ prim ) 4 . (21b) 

Bounds (20ab) would certainly be less compelling if the Fly's Eye event of Fig. 4 were 
unique. But, luckily, this is not the case, as the Pierre Auger Observatory [29] has already 
seen an event (ID No. 737165) with a similar energy of 2 x 10 11 GeV [30]. 

Let us end this subsection with a somewhat peripheral remark. The extremely small 
numbers of order 10~ 38 msik/(2x 10 22 GeV) appearing in the Cherenkov-like bound 
(20b), which traces back to Eq. (1.14) of Ref. [27], appear to be incompatible with a 
quartic-term mass scale Mqg2 ~6x 10 10 GeV claimed [31] to be a possible explanation 
of time-dispersion effects observed in the July 9, 2005 gamma-ray flare from Mkn 501 
by the MAGIC telescope (see, in particular, Fig. 7 of Ref. [32]). Even though we have 
doubts as to the validity of this possible detection of "quantum-gravity" effects, 2 it is 
instructive to see how, in principle, astrophysical data could give more than just bounds 
on the possible small-scale structure of space. 



2 Note also that the June 30, 2005 flare, for the two medium-energy bands 0.25 — 0.60 TeV and 0.60 — 
1 .20 TeV shown in Fig. 6 of Ref. [32], appears to have a time-dispersion behavior opposite to that of the 
July 9, 2005 flare, as shown in Fig. 7 of the same reference. 
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3.2. UHECR Bounds - Model 2 



Next, obtain Cherenkov bounds [2, 3] for the nine parameters a of Model 2. Start 
with the 15 Auger events [30] of Table 1 on the next page, where the event identifica- 
tion number is given in column 1, the energy E pr - m of a hadron primary in column 2, the 
shower- maximum atmospheric depth X mSLX in column 3, and pseudo-random event direc- 
tions with right ascension RA' G [0,360°] and declination DEC' G [-70°, 25°] in column 
4. 3 These are high-quality events, having been observed in the "hybrid" mode (Fig. 5) 
of the Pierre Auger Observatory [29]. Specifically, the X max values give information on 
the primary particle type, as discussed in Ref. [3] and references therein. 

As mentioned in Sec. 2.2, vacuum Cherenkov bounds can only be obtained for 
parameters a in domain (18), for which the phase velocity of light is less than the 
maximal attainable velocity c of the charged particle. Now determine a hypersphere S 8 
of radius a in this subspace, so that for each a G S 8 the Cherenkov threshold condition 
(19), with expression (14) inserted, is violated for at least one event from Table 1. 

The excluded domain of parameter space then corresponds to the region on or outside 
this hypersphere S~. Namely, for a positive integrand of the ramp function on the right- 
hand side of (14), the first term on the right-hand side is multiplied by a factor 1/A 
under scaling of a — > X a (and e — > X e) with X > 1, so that inequality (19) is violated 
for X > 1 if it is already for X = 1 . The events from Table 1 give at the 2c level: 

deluded = »cauL"| ,n \&\>»}, (22a) 

■ 53 3xl0 " 8 (l6§v7?) 2 ' < 22b > 

where \a\ is the standard Euclidean norm (12) and where the mass of the primary 
charged particle has conservatively been taken equal to that of oxygen (most primaries 
in the sample are expected to be protons and helium nuclei). 

The resulting UHECR Cherenkov bounds for nonbirefringent modified-Maxwell the- 
ory can be described as follows: considering deformation parameters a , . . . , a 8 with 
a corresponding phase velocity of light less than the maximal attainable velocity of 
charged particles, each of these nine parameters must separately have a modulus less 
than the value a given by (22b). Improving slightly by considering all nine parameters 
simultaneously, the new 2o bound from astrophysics is given by 

|a| 2 = £(a') 2 <(3xl0- 18 ) 2 , (23) 
i=o v 7 

with MpHm simply set to 16 GeV/c 2 . 



The directions of these 15 events have not yet been released by the Pierre Auger Collaboration and, for 
this reason, fictional directions have been given in Table 1 (with primes alerting to their nonreality). The 
bounds of this subsection are based on the numbers given in Table 1, but we are confident that the same 
bounds are obtained from the actual event directions when they are made available by the Pierre Auger 
Collaboration. For the bounds obtained here, it only matters that there are 15 UHECR events with more or 
less equal energies E prim w 15 EeV and more or less random directions over a significant part of the sky, 
the precise association of energy and direction being irrelevant. 
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x direction cosine 

FIGURE 4. Pointing directions of the 22 photomultiplier tubes which triggered in connection with the 
Fly's Eye event [26] of October 15, 1991 at 7:34:16 UT. The pointing directions are shown projected 
into the x-z plane, where the jt-axis points east, the y-axis north, and the z-axis upward. The triggered 
phototubes have positive y-components. The dashed line indicates the plane defined by the shower axis 
and the detector. See Ref. [26] for further details (figure reproduced by permission of the AAS). 



TABLE 1. "Hybrid" events recorded by the Pierre Auger Obser- 
vatory [29, 30]; see Sec. 3.2 in the main text for further details. 



ID No. 


E prira [EeV] 


Xmax [gem 2 ] 


(RA', DEC') [deg] 


668949' 


18 


765 


(356 , -29) 


673409' 


13 


760 


(344 , -62) 


828057' 


14 


805 


(086 , -34) 


986990' 


16 


810 


(152,-33) 


1109855' 


17 


819 


(280 , -30) 


1171225' 


16 


786 


(309 , -70) 


1175036' 


18 


780 


(228, +17) 


1421093' 


27 


831 


(079, +13) 


1535139' 


16 


768 


(006 , -62) 


1539432' 


13 


787 


(153 ,-15) 


1671524' 


14 


806 


(028 , -63) 


1683620' 


21 


824 


(024 , -23) 


1687849' 


17 


780 


(031 ,-23) 


2035613' 


12 


802 


(079 , -08) 


2036381' 


29 


782 


(158,-03) 




FIGURE 5. The hybrid nature of the Pierre Auger Observatory [29] provides for two independent ways 
to study extensive air showers, namely, by an array of surface water-Cherenkov detectors and a collection 
of air-fluorescence telescopes [see http : / /www. auger . org/]. 
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Three remarks are in order. First, bound (23) holds for the reference frame in which 
the cosmic-ray energies are measured. This frame is essentially given by the rest system 
of the Sun. Second, the restriction of a to domain (18) makes the Cherenkov bound 
(23) effectively "one-sided," as discussed in the third paragraph of Sec. 3.1. To specify 
the precise domain in the 9-dimensional parameter space requires, however, some care. 
For example, in the 2-dimensional subspace with only Or = a cos0 and a 1 = a sin0 
nonzero, the allowed domain from (23) is given by an open disk segment with —n/4 < 
(j) < it I A and < a < 3 x 10~ 18 . Third, consider the restricted isotropic model with 
only parameter cc° nonzero, which was already discussed in the paragraph containing 
the exact tree-level results (15)— (17). Then, the single Fly's Eye event [26] with E pr \ m as 
300 EeV (or the single 200 EeV Auger event [30] mentioned in Sec. 3.1) gives a one- 
sided bound on a at the 10~ 23 level; see Eq. (C6) of Ref. [2]. This single a bound 
corresponds, in fact, to the negative part of the previous bound (20a). 

At this moment, it may be of interest to recall the best bounds from laboratory 
experiments [33, 34, 35]: 

• for the spatially isotropic nonbirefringent deformation parameter oc°, there is a di- 
rect bound at the 10~ 7 level [33(c)] and an indirect (electron anomalous-magnetic- 
moment) bound at the 10~ 8 level [34(c)]; 

• for the eight other nonbirefringent parameters Of 1 , . . . , a 8 , there are direct bounds 
at the 10~ 12 to 10~ 16 levels [35]. 

Observe that, for Earth-based cavity experiments with an orbital velocity v© around 
the Sun, the parity-odd parameters a 1 , for / = 1,2,3, typically enter the measured 
relative frequency shift with a factor (v©/c) 10~ 4 , which reduces the sensitivity for 
these three parameters. In these and other experiments, the parameter 05° even enters 
with a quadratic boost factor (precisely which velocity plays a decisive role depends 
on the experimental setup) and the sensitivity is reduced also for this parameter. See 
Refs. [15, 33, 35] for further details. 

Astrophysics and laboratory bounds on a are not directly comparable, the latter 
having all the benefits of repeatability and control. But the astrophysics bounds at the 
10~ 18 level are certainly indicative. Moreover, the UHECR Cherenkov bounds on a can 
be expected to drop to the 10~ 23 level in the coming years [3]. 



3.3. UHECR Bounds - Implications 

In Sec. 3.1, we have obtained Cherenkov bounds for the effective defect excluded- 
volume factor F and the effective defect size b in the modified photon dispersion relation 
(3b): 

< F < 10~ 23 , (24a) 
< b < 10~ 26 m « hc/(2x 10 10 GeV), (24b) 

where the last bound follows from (20b) with 04 F = 3 x 10~ 23 inserted. But even the 
"weakened" bound (24b) is already quite remarkable compared to the previous bound (1) 
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mentioned in the Introduction and, more generally, to what can be achieved with particle 
accelerators on Earth (the upcoming Large Hadron Collider at CERN will have a proton 
beam energy E p ^l x 10 3 GeV). In fact, bound (24b) appears to rule out so-called TeV- 
gravity models; see Endnote [49] of Ref. [1], which contains further references. 

Still, the most interesting bound may very well be (24a), as it rules out a single- 
scale classical spacetime foam having approximately equal effective defect size b and 
separation /, with an excluded- volume factor F = (b/l) 4 ~ 1. Moreover, this conclusion 
holds for arbitrarily small values of the effective defect size b, as long as a classical 
spacetime makes sense. 

More generally, also the nineteen Lorentz-violating deformation parameters of the 
modified-Maxwell theory (6) are strongly bounded, as discussed in Sec. 3.2: 

< \K^ vpa \ < 1(T 18 , (25) 

where, for the sake of argument, nine "one-sided" Cherenkov bounds from (23) have 
been made "two-sided" (the even stronger astrophysics bounds [15] on the ten bire- 
fringent parameters are already two-sided). A similar bound has been obtained for the 
stochastic model mentioned a few lines below (1 1); see, in particular, bound (1 lb) from 
Ref. [17(b)]. 

A priori, one expects 0(1) effects for (24a) and (25) if there is some kind of nontrivial 
small-scale structure of classical spacetime. This expectation is based on an early obser- 
vation by Veltman [36] that, if a symmetry (for him, gauge invariance) of the quantum 
field theory considered is violated by the high-energy cutoff A (or by a more funda- 
mental theory), then, without fine tuning, the low-energy effective theory may contain 
symmetry-violating terms which are not suppressed by inverse powers of the cutoff en- 
ergy A. The same observation holds for Lorentz invariance, as emphasized recently by 
the authors of Ref. [37]. 

Contrary to these expectations, the experimental bounds (24a) and (25) seem to indi- 
cate that Lorentz invariance remains perfectly valid down towards smaller and smaller 
distances. In fact, this conclusion would appear to hold down to distances at which the 
classical-quantum transition of spacetime occurs (generally considered to be equal to 
the so-called Planck length). 

In the next section, we will try to make a first step towards understanding the ro- 
bustness of Lorentz invariance, starting from a physical picture based on a dynamic 
spacetime and quantum uncertainty (entirely in the spirit of Einstein and Heisenberg). It 
may be of interest to mention one other possible explanation in the context of so-called 
"emerging-gravity" models [38, 39, 40, 41], which relies on the existence of a trans- 
Planckian energy scale in a fundamental Lorentz-violating theory with a topologically 
protected Fermi point [42]. But this avenue will not be pursued further here. 

4. QUANTUM-SPACETIME CONJECTURES 

The main conclusion from Sec. 3.3 is that Lorentz invariance appears to remain valid 
down towards smaller and smaller distances, even down to distances at which the 
classical-quantum transition of spacetime is believed to occur. That distance is usually 
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taken to be the Planck length [7, 43, 44], 



Planck = yhG/c 3 « 1.6 x 1CT 35 m « »c/(l.2 x 10 19 GeV) . (26) 

The question, however, is whether or not quantum spacetime effects show up only at 
distances of the order of the Planck length. Perhaps a quantum spacetime foam [6, 7, 
8, 9] could arise primarily from gravitational self-interactions which need not involve 
Newton's constant G describing the gravitational coupling of matter (similar to the case 
of a gas of instantons in Yang-Mills theory). 
This brings us to the following suggestion: 

Conjecture 1 Quantum spacetime has a fundamental length scale I which is essentially 
different from the length /pi an ck as defined by (26). 

Concretely, this length / could have a different dependence on G than the Planck length. 
In the next two subsections, we will argue that the standard theory appears to leave some 
room for a new fundamental constant / and, in the third subsection, we will present a 
further conjecture as to what physics may be involved. 

4.1. Generalized Action 

Our starting point is Feynman's insight [45] that the quantum world of probability 
amplitudes for matter is governed by the complex phase factor exp(i J^m) with a real 
phase J^m = Sy\/h expressed in terms of the classical matter action 5m and the reduced 
Planck constant h = h/2n. 

As the merging of quantum mechanics and gravitation is far from understood, it may 
be sensible to consider a generalized dimensionless action for "quantum gravity" (or, 
better, "quantum spacetime") as probed by classical matter: 

^ grav = I d4 *v^M (*(*) + 2X ) + / d 4 *v^«T^M lass W , (27) 

with the Ricci curvature scalar R(x), a new fundamental length scale Z, and a nonnegative 
cosmological constant A > 0. 

Two remarks may be helpful. First, expression (27) is to be used only for a rough 
description of possible quantum effects of spacetime, not for those of the matter fields 
which are considered to act as a classical source with coupling constant G. Second, the 
special case I 2 = /p| anck reduces the dimensionless action (27) to the standard expression 
[9], with the Einstein-Hilbert integral multiplied by the factor —c 3 /(16kG) and the 
matter integral by l/h. Here, however, we wish to explore the possibility that / is an 
entirely new length scale independent of /pianck (Conjecture 1). 

The overall factor l~ 2 in (27) is irrelevant for obtaining the classical field equations, 

R» v (x) - \ g^{x)R{x) - ?Lg» v (x) = -%KGT^{x) , (28) 

where the matter energy-momentum tensor 7^ (x) is defined by an appropriate func- 
tional derivative of / d 4 x y/\g\ ^M lass - But, following Feynman, the physics of a genuine 
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quantum spacetime would be governed by the complex phase factor exp(i^ grav ) and 
the overall factor of l~ 2 in (27) would be physically relevant. Just to be clear, we do not 
claim that a future theory of quantum gravity must necessarily be formulated as a path 
integral but only wish to use the Feynman phase factor as a heuristic device in order to 
say something about the typical length scale involved. 

The expression (27) for the generalized quantum phase also suggests that, as far as 
spacetime is concerned, the role of Planck's constant h would be replaced by the squared 
length I 2 , which might loosely be called the "quantum of area" (see, e.g., Ref. [46] for a 
calculation in the context of loop quantum gravity and Refs. [47, 48] for general remarks 
on a possible minimum length). Planck's constant h would continue to play a role in the 
description of the matter quantum fields. 

However, with h and I 2 being logically independent, it is possible to consider the 
"limit" h — > (matter behaving classically) while keeping I 2 fixed (spacetime behaving 
nonclassically). Even if the numerical values of the lengths / and /pi an ck will turn out to 
be close (or equal) in the end, it may be conceptually interesting to consider phases of 
the theory with ratios ///pi anc k very different from unity. 

4.2. Fundamental Constants 

The fundamental dimensionful constants encountered in the previous subsection are 
summarized in Table 2 on the next page. In this contribution, we mainly consider the 
second and third columns of Table 2 (removing the parentheses around G) and leave 
the rigorous treatment of all columns to a future theory. In this subsection, however, we 
already make some speculative remarks on a possible new theory encompassing all of 
the physics arenas appearing in Table 2. 

In that future theory, "classical gravitation" may perhaps be induced [49] by combined 
quantum effects of matter and spacetime, giving Newton's gravitational constant 

G=Cc 3 l 2 /h, (29) 

with a calculable numerical coefficient £ > 0. 4 In this way, the "large" classical coupling 
constant G would be the ratio of "small" quantum constants I 2 and h. Practically, the 
numerical coefficient £ gives the ratio lp\ anc ^/l 2 - 

The dimensionful constants h, c, and I 2 of Table 2 correspond to, respectively, the 
quantum of action, the limiting speed of a massive particle, and the quantum of area 
(up to an as yet undetermined numerical factor). In addition, each of these constants 
is believed to play "the role of a conversion factor that is essential to implementing a 
profound physical concept" (quoting from Ref. [44]). 

The role of the individual constants c and h is clear: the first converts time into length 
according to relativity theory (cT = L, where physical quantities to be converted are 
denoted by capital letters) and the second converts frequency or inverse time into energy 



There may or may not be a connection with the "emergent-gravity" models mentioned in the last 
paragraph of Sec. 3.3. See also Ref. [50] for an interesting comparison to hydrodynamics. 
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TABLE 2. Fundamental dimensionful constants of nature, in- 
cluding the hypothetical quantum of area I 2 . The gravitational con- 
stant G appears between brackets, as it can, in principle, be ex- 
pressed in terms of the other constants; see (29) in the main text. 

quantum matter classical relativity quantum spacetime 

ft c (G) I 2 



according to quantum theory (hQ. = hln/T = E). From Einstein's theory of gravity 
and relation (29), it would seem that the combination I 2 /(he) acts as the conversion 
factor between space curvature (with dimensions of inverse length square) and energy 
density (with dimensions of energy over length cube), so that 1/L 2 ~ £ I 2 /(he) E /L? , 
with the numerical coefficient £ from (29). Replacing the energy E by Me 2 , these 
conversion factors are summarized in Fig. 6. Starting with LENGTH on the left of this 
figure and going around clockwise, the required number of fundamental constants is 
seen to accumulate steadily. 

In the previous paragraph, we have found the conversion factor between mass and 
length by relying on Einstein's theory of gravity. Incidentally, the simplest way to get 
this conversion factor may be to use the standard result for the Schwarzschild radius 
of a central mass M, i?Schw = 2GM/c 2 , and to replace G therein by (29). For the 
corresponding arrow in Fig. 6, this would imply that the conversion factor l 2 c/h comes 
together with the pure number £ appearing in (29). However, it might be interesting, as 
discussed in the last paragraph of Sec. 4.1, to consider a possible phase of the theory 
with C = and still keep the linear relation between mass and length. 

For £ = 0, the physical interpretation of the conversion factor l 2 c/h in Fig. 6 is 
entirely open (hence, the dotted arrow in the figure). A trivial suggestion (definitely not 
yet a "profound physical concept") would run as follows. Given the fundamental length 
scale / and the constants c and h, a corresponding energy density is defined as p = I? c/7 4 . 
Then, one possible physical interpretation would be that the mass M of an "elementary 
particle" is associated to a characteristic length L in such a way that the fundamental 
energy density p over a volume V obtained by multiplying L with the minimum area I 2 , 



M[\/c] 

LENGTH " 




M[l 2 c/h] 



FIGURE 6. Fundamental dimensionful constants h, c, and I 2 from Table 2 in the role of conversion 
factors, with numerical coefficients omitted. The operation M[z] stands for multiplication by z and the 
operation / for inversion. Hence, time is first inverted (to give a frequency) and then multiplied by h/c 2 
(to give a mass). The linear relation between mass and length is discussed in Sec. 4.2 of the main text. 
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V = LI 2 , returns the original mass, M = Vp. In this way, the associated length L of a 
mass M would correspond to the largest linear spread a fiducial volume V =M/p could 
have, provided the area of any cross section of the volume V is not allowed to drop under 
the value I 2 . The heuristic idea behind this length L is that an inserted material mass M 
distorts the intrinsic foam-like structure of space. 

Observe that the mass dependence of the associated length L = Ml 2 c/fi is opposite 
to that of the standard Compton (reduced) wavelength %q = Ti/(Mc). For the known 
elementary particles, with M of the order of 10 2 GeV/c 2 or less, and I below the value 
quoted in (24b), the mass-induced fuzziness scale L of a particle will be very much 
smaller than the corresponding Compton wavelength Xq. Specifically, the ratio of these 
two length scales is given by the inverse ratio of the corresponding squared energies, 

L/Xc = {Mc 2 /(hc/l)) 2 . 

At this moment, it may be helpful to list all the length scales possibly relevant to 
an elementary particle with mass M (in order of increasing size, taking £ < 1/2 and 
temporarily setting c = h = 1): the Schwarzschild radius 7?s c hw = 2£/ 2 M, the mass- 
induced fuzziness scale L = l 2 M, the fundamental quantum-spacetime length scale 
/, and the Compton wavelength Xq = \/M. In principle, the mass-induced fuzziness 
scale L (or perhaps even the larger intrinsic fuzziness scale /) could lead to form-factor 
effects in high-energy scattering processes and possibly affect the renormalization of the 
standard-model quantum field theory. A minimum length may even play a crucial role 
in obtaining the usual rules of quantum mechanics, as discussed in Ref. [51]. 



4.3. Vacuum Energy Density 

Let us return to the generalized action (27), with Newton's gravitational constant G 
appearing explicitly but without Planck's constant fi. If Conjecture 1 holds true, the 
question arises as to what type of physics determines the length scale /. Here, we assume 
this physics to be independent of h, otherwise the discussion would be more along the 
lines sketched in the previous subsection. One possible answer would then be given by 
the following suggestion: 

Conjecture 2 The quantum-spacetime length scale I is related to a nonvanishing cos- 
mological constant or vacuum energy density. 

For the case of the early universe, with a vacuum energy density p vac = E* ac / (he) 3 
in the matter part of the action (27) and X = in the geometric part, it can be argued [4] 
that the following approximate relation holds (h appears only as an auxiliary constant): 

/ ~ c 2 /^Gp— = hcE M /El c ^ l(T 29 m ( (^f^ ( 3 °) 

where the Planck energy scale is given by fipianck = ^ c /^Planck an d the numerical value 
for E yac has been identified with the "grand-unification" scale suggested by elementary 
particle physics [5, 52]. As anticipated in the sentence under Conjecture 1, the length 
scale / as given by the first mathematical expression on the right-hand side of (30) has a 
G dependence which is different from that of the standard Planck length (26). 
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If (30) holds true with /pi anc kA ~ 10 6 ' it i s perhaps possible to have sufficiently rare 
defects left over from the crystallization process of an initial quantum spacetime foam to 
a classical spacetime. 5 With the effective defect size b set by /pi an ck (matter related) and 
the effective defect separation / set by / (vacuum related), these spacetime defects would 
give the following excluded-volume factor in the modified photon dispersion relation 
(3b): 

F=(b/l) 4 I 1(T 24 , (31) 

which is close to saturating the current UHECR bound (20a). As mentioned at the 
end of Sec. 3.1, a possible detection of Lorentz- violating effects (having ruled out 
conventional explanations) would be an entirely different matter than setting better and 
better bounds on the small-scale structure of spacetime, even though the latter type of 
"null experiments" can also be of great importance (think of the role of the Michelson- 
Morley experiment for the discovery of special relativity [54]). 

5. SUMMARY 

The main phenomenological conclusion of the first part of this review (Sees. 2 and 3) 
is that quantum spacetime foam, if at all real, appears to have given rise to a classical 
spacetime manifold which is remarkably smooth. Using astrophysics data, this result 
can be quantified as follows: the defect excluded-volume factor in the modified proton 
and photon dispersion relations (3ab) is bounded by F < 10~ 23 <C 1 and the Lorentz- 
violating parameters in the modified-Maxwell theory (6) by \ K^ vpa \ < 10~ 18 <C 1. 

Prompted by this phenomenological conclusion, a theoretical suggestion has been 
advanced in the second part (Sec. 4), namely, that the quantum theory of spacetime 
may have a fundamental length scale / conceptually different from the Planck length 
(Conjecture 1) and possibly related to vacuum energy density (Conjecture 2). With two 
length scales present, I and /pi an ck> it is perhaps possible to satisfy the tight experimental 
constraint on F. 

Still, it is safe to say that the quantum origin of classical spacetime remains a mystery. 
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5 It is likely that a proper understanding of this crystallization process requires new ideas about the 
quantum measurement problem and possibly even a mechanism for objective state reduction [53]. 
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